Abstract 

This paper continues the study of the properties of an accretion disc rotating around a non-baryonic 
(assumed super-massive) compact object. This kind of objects, generically known as boson stars, were 
earlier proposed as a possible alternative scenario to the existence of super-masive black holes in the 
center of every galaxy. A dilute boson star has also been proposed as a large part of the non-baryonic 
dark matter, flattening galactic rotational velocities curves. In this contribution, we compute the profile 
of the emission lines of Iron; its shape has been for long known as a useful diagnosis of the space-time 
geometry. We compare with the case of a Schwarzschild black hole, concluding that the differences are 
observationally distinguishable. 
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1 Introduction 

a: 

^ It has long been believed that accretion onto massive black holes is the energy source of active galactic 
^ nuclei (AGN). Nearby galaxies, as the remnants of AGNs, should then harbor massive black holes in their 
• t-h centers [||. Recent observations have indeed revealed the existence of massive dark objects, believed to be 
^ chose massive black holes, in most, if not all, nearby galactic nuclei [§, [|. However, those observations 
only reach a physical radius 10000 times larger than that of the event horizon of the putative black hole. 
As pointed out in Ref. [Q], from a dynamical point of view, a single massive boson star is also compatible 
with all measurements of star kinematics around the massive dark object in the most studied center, that 
of Milky Way, without requiring the massive dark object to be a black hole. A proof of the existence of 
massive black holes requires unambiguous evidence of the event horizon. It is then important to address 
the issue of the emissivity properties of matter rotating around different central objects. Ref. || made 
a first step in this sense; to that paper we refer the reader for further discussions and a more complete 
list of references on boson star physics and related problems. The emissivity properties are much more 
important than the kinematic studies, since they can provide a more direct way to actually distinguish 
between different central objects ||. 

X-ray studies have revealed the most prominent spectral feature coming from the central regions of 
some active galaxies, an enormous Doppler and gravitational-shifted broad Fe Ka line, which, for the first 
time, showed the existence of a strong gravity regime [0, (8|. This has been taken as the strongest evidence 
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for the existence of super-massive black holes, or even as a sensitive diagnosis of the black hole spin. But, 
how would the profile of the emission lines be for those produced in discs rotating around putative boson 
stars? Would them present sufficiently different shapes as to make them distinguishable in forthcoming 
observations? These are the questions we address in this paper. 

2 Photon propagation in the metric of a boson star 
2.1 Space-time 

The space-time generated by a stationary and spherically symmetric object can be described by the metric: 

ds 2 = B(r)dt 2 - A(r)dr 2 - r 2 (dQ 2 + sin 2 6># 2 ), (1) 

where A(r) can be written as 1/(1 — 2GM(r)/rc 2 ), c is the speed of light, G is the gravitation constant 
and M[r) is the mass enclosed within a radius r. The simplest solution for Eq. (jl]) is that representing a 
Schwarzschild black hole, for which B(r) = 1 — 2GM/rc 2 , and M(r) is a constant equal to the mass of the 
black hole. 

Other stellar-like systems described by Eq. (|l|) are those known as boson stars. Details of boson star 
solutions can be found in the reviews quoted in Ref. ||]. For a boson star, A(r) and B(r) can be calculated 
only numerically (see 0] for explanations, the program herein used is an adaptation of that mentioned in 
Refs. 0). 

The total mass, the radius and all other physical parameters of the boson star are scaled by the mass 
of the constituent boson, m. For numerical purposes, then, we can work with dimensionless units. The 
radius in a dimensionless form will be x = mr, and the mass, as well, will be transformed into a dimen- 
sionless quantity, M = M(x)Mp Unck /m, where Mpi anc k = \Jfic/G, h is the Planck constant. This scheme 
was already proposed in the first papers on boson stars p| , |l2| , |l~3|| . For example, if the massive dark 
object with mass 2.6 x 1O 6 M0 in the Galactic Center is a boson star, then the constituent boson mass is 
m = 2.81 x 10~ 26 GeV. The distance from the center in unit of parsec is r[pc] = x / m[GeV] 6.38 x 10~ 33 , 
the mass in millions of solar mass is M[10 6 ] = M(x)/m[GeV] 1.33 x 1(T 25 , and the gravitational radius 
r g = GM/c 2 is written as x g = 0.53256 in dimensionless units. The boson star structure is numerically 
computed to obtain M(x) and B{x) Q. Two notes are, however, in order. Firstly, these numbers rely in 
that we are using a non self-interacting boson star, for simplicity and because the behavior is generic, but 
should be consistently changed if this is not the case Q. By going towards a self- interacting model, the 
mass of the constituent boson needed to produce the same overall mass increases. Secondly, the use of the 
mass of the galactic center region is done just to provide a numerical example, we are not claiming the 
possible existence of a boson star in the center of the Galaxy, since apparently the emissivity properties 
would conflict with observations (see Ref. ||), or at least we know that that is the case for the simplest 
accretion disc models. 

To conveniently compare the metric generated by a non-rotating boson star with the Schwarzschild 
black hole metric, we re-scale the total mass of the boson star to be M = 1. We do so by just dividing 
x/0. 53256 and replacing Mix) by M(x)/0. 53256. Figure |l| shows B(r) as a function of r for both boson 
star and Schwarzschild black hole. B(r) for both cases are almost exactly the same when r > 15r g , while 
it slowly decreases to a constant at the center of a boson star rather than rapidly decrease to zero at the 
horizon (r = 2r g ) of a Schwarzschild black hole. To a distant observer, thus, a photon emitted from the 
inner region of a disc extending within a boson star (r < 15r g ) suffers much less gravitational redshift than 
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Figure 1: The metric parameter B(r) as a function a radius (in unit of the gravitation radius r g = GM/c 2 ). 
The solid line represents a boson star and the dash line represents a Schwarzschild black hole of equal mass. 
The main difference appears in the inner region, with r < 15r g . 



the one emitted from the disc surrounding a Schwarzschild black hole. 



To simplify the following calculations, we adopt a simple Logistic function to fit both M(x) and B(x) 
for x < 50, 

M(x)[ovB(x)] = f 1 ~^ p +A 2 . (2) 

For the mass Mix), the coefficients are A x = 0.01039, A 2 = 0.53533, x = 5.49809 and p = 3.6201; while 
for B(x) the coefficients are A 1 = 0.78236, A 2 = 0.99793, x = 8.17558 and p = 2.09499. Both fits are 
excellent (x 2 < 10 -5 ) and valid throughout all space-time (see Fig|2|); they much improve those presented 
in Ref. For x > 50, M{x) and B{x) are almost exactly the same as those obtained in the Schwarzschild 
black hole case, i.e. M(x) = 0.53256 and B(x) = 1 - 2M(x)/x. 



2.2 Orbits 



The Lagrangian of the geodesies in the space-time with a line element shown by Eq. (|T|) can be written as 
(cf. |3, pages 342-347 — ||, pages 897-902): 



2C = B(r)t 2 - A(r)r 



>',2 



r 2 (9 2 + sin 2 i 



(3) 



where the dot denotes differentiation with respect to an affine parameter A along the geodesies. Together 
with the assumption of stationary and axisymmetric space-time, the Euler-Lagrange equations give the 
simplified photon momentum equations as (see analogous derivations of the momentum equations for black 



hole metric in Ref. [TL4L 15], pages as above): 



dt 
dX 



B(r) s/R' 



(4) 
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Figure 2: Logistic function fits to both, B(x) and the mass. The fit is actually superposed to the numerical 
curves in both cases. 
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where the functions and R are respectively given by 

= r]-i 2 cot 2 e 
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A{r)B{r) A(r) 
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(5) 
(6) 
(7) 

(8) 
(9) 



and where r/ = K/E 2 and £ = L/E are constants of motion related to the angular momentum L, the energy 
at infinity E, and the Carter's constant K. = pg + L 2 cot 2 9. p$ is the ^-component of photon momentum 



(cf. the photon momentum equations in the space-time of a Kerr black hole [14, 15], pages as above). 

The path of a particular photon is completely described by the motion constants r] and £, which are in 
fact related to the polar and azimuthal angles a and (5 on the sky as seen by a distant, locally non-rotating, 
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observer who receives the photon. The motion constants r\ and £ are then given by Q 

/ r sin 9 sin a sin j3 \ 

v 7W) L„ ' 



B(r) 



sin 2 a (l — sin 2 sin 2 f3 



The equation governing the projection of a photon orbit in (r, #) plane is 



dr 
1 



71' 



(10) 

(ii) 

(12) 



where r, and 9% are the initial values of r and 0. The ^-integration, i.e. the right hand of Eq. (O), gives 

dO 
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To enable the integration of a given photon away from the observer to the disc surface, it is necessary 
to determine the initial sign of dr/dX and d9/dX which are in fact given by (see Fig. la in Ref. Jl6|): 



(16) 
(17) 
(18) 
(19) 



Armed with Eqs. (Q) to (|T9|), we can trace the trajectories of photons that propagate from the observer 
at (r , 9 , <p ) back to the initial emitting position (r c , vr/2, <j) e ) in the local rest frame of matter moving 
in the plane of the disc. 

1 See the derivation in Ref. Jl6|, 01 for the metric of Kerr black hole. The motion constants could be obtained by replacing 
the Kerr metric in Karas et al. with the much simple boson star metric. See Fig. la in Ref. jl6| for a geometric 

interpretation of the angles involved. 
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The sign of dr/dX (or d9/dX) changes at the turning point where R = (or = 0). 
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2.3 Disc 



In analogy to the standard geometrically thin and optically thick disc around a black hole, such simplified 
accretion disc now surrounding a massive boson star was first discussed in Ref. ||. It is assumed that 
the disc material is in perfect Keplerian motion which gives the following components for the 4-velocity 
(analogous to the derivation for equatorial disk around Kerr black holes in Ref. 

(20) 




(21) 

(22) 
(23) 

where B' (r) = dB(r)/dr. The rotating velocity of the disc material measured by a static observer is given 
by 

y(*)=rd B ' [r) , ■ (24) 
V 2B(r)-rB'(r) 1 ' 

The coordinate-velocity in Eq. (21) is related to angular velocity by ur = -u'fl^ Note that in the potential 
of a non-rotating boson star, there are circular orbits for every possible value of the radial coordinate, 
including those which are inside the boson star ||. This is clearly different from the black hole case, where 
the innermost stable circular orbit only extends down to r = 6r g , in the case of a Schwarzschild black hole. 
As seen from Figure ||, at large radii (r > 25r 9 ), the rotating velocity of the disc surrounding a boson star is 
the same of that disc surrounding a Schwarzschild black hole. At the innermost part of the disc (r < 25r g ), 
however, differences are not negligible. In the case of a boson star, with decreasing radius, the rotating 
velocity of the disc increases to a maximum value and then decreases to zero at the star center. In the 
case of a Schwarzschild black hole, the rotating velocity of the disc increases monotonically until the radius 
approaches the last stable orbit (r = 6r g ). To the view of distant observers, then, this difference results 
in less Doppler-shifted photons observed in the putative case of a boson star-accretion disc (BS-AD), as 
compared to the Schwarzschild black hole accretion disc (SchBH-AD). 

2.4 Numerical implementation 

We shall assume that 1) the disc material is not interacting with the bosons, 2) the distant observer views 
the disc at some inclination angle o (with O = 0° corresponding to a face-on disc). Operationally, we 
define the distant observer to be located at r Q = 10 4 r g and fa = 0°. We then integrate the four photon 
momentum equations from the observer to the accretion disc using motion constants £ and r/. The redshift 
factor of a photon propagating from the disc to the observer is given by 

U Q Po • Up , , 

9 = — = , (25) 

v e p c • u e 

where v Q and v e are the frequencies of the photon received by the distant observer and the photon emitted 
in the rest frame of disc material, respectively; p Q is the photon momentum received by the distant observer 
and p e is the photon momentum in the rest frame of the emitting material. The disc region between r out 
and r; n is divided into 1000 x 1000 elements (ji — ► r\ + dr, fa — > fa + dfa). The solid angle subtended at 
the observer's sky by each element is also calculated by counting the photons impinging onto the element 



2 Note that there is a typing error in the first of the Eqs. (3) of Ref. H: x should be in the denominator. 
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Figure 3: The rotating velocity of the disc material in the view of a static observer as a function a radius 
(in unit of the gravitation radius r g = GM/c 2 . The solid line represents a boson star and the dashed line 
represents a Schwarzschild black hole. 



Similar to the case of an accretion disc surrounding a black hole, the disc rotating around a boson star 
may also be illuminated by some nearby hard X-ray source due to similar mechanism, for example, inverse 
Compton scattering of soft thermal photons by high energy electrons in the corona above the disc. The flu- 
orescent Fe Ka line photons driven by this illumination may directly propagate to distant observers without 
absorption. We assume, as it is usual, that the X-ray illumination is axisymmetric and the emissivity of 
Fe Ka line only varies with radius as a power law e(r) oc r~ a in the rest frame of emitting material, where 
a is the emissivity index. We already know that the viscosity dissipated power per unit area approaches a 
maximum value at the radius of ~ 5r g and then rapidly decreases as the accreted material goes into the 
star center ||. This is quite different from the case of black hole accretion disc, where the radiation from 
the region within the marginally stable orbit (in projection) is insignificant. 

In order to see the difference between the profiles of the emission lines coming from a BS-AD and from 
a SchBH-AD, respectively, we consider two cases. In the first one, for a BS-AD, the line is assumed to be 
emitted from an annulus between r; n = 6r g and r ou t = 30r g . In the second case, the line will be assumed 
to come from an annulus between r; n = 2r g and r ou t = 30r g . At all times, for a SchBH-AD, the line 
is assumed to be emitted from annulus between rj n = 6r g and r out = 30r g . Physical considerations can 
be put forward to justify these values of inner and outer limits of the disc. In the case of BS-AD, r- m is 
not chosen to be closer to the center (or 0) because we are adopting a power-law emissivity law which 
diverges at the star center. We note, however, that the line emission within r < 2r g would be negligible if 
the emissivity law is similar to the viscosity dissipated power as a function of radius. We set r ia = 6r g in 
the case of SchBH-AD since the disc only extends to the marginal stable orbit. The line emission from the 
disc material within the marginal stable orbit is generally assumed negligible, and if this line flux is even 



(we refer the reader to Refs. [17, [H| for further details). 
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observable, as Reynolds and Begelman (20] have pointed out, it would generally emerge in the extremely 
red wing of the line profile, what would make the difference between the profiles of Fe Ka generated from 
the BS-AD and from the SchBH-AD that we present in the next Section larger. 



The observed line flux distribution is dF(E Q ^ s ) = I(E b s )dQ, where d£l is the solid angle subtended by 
the disc in the observer's sky and I(i5 bs) is the specific intensity received by the observer. The specific 
intensity I(E em ) of the line emitting in the rest frame of disk material is assumed to be isotropic and 
monochromatic I(E em ) = e(r)5(E cm — Eq), where Eq (= 6.4 keV for iron Ka line) is the energy of the 
line photons in the rest frame of the emitting material and gE em = E Q ^ S . From the relativistic invariance 
I u /v 3 (see explanations in Ref. [15], pages 897-902), we have /(-E bs) = {v Q / v e ) 3 1 (E em ) = g 3 I(E em ) = 
g 4 e(r)<5(-Eobs — qEq), and thus the total observed flux distribution of the line is given by 

/•r-out 

F(E ohs ) oc / e(r)<7 4 5(£ obs - gE )dn oha , (r in < r < r out ). (26) 



This is what is plotted in the graphs of next Section. 



3 Line Profiles 

Theoretical emission line profiles produced in the inner region of accretion discs around black holes have 
been calculated by several authors, see for instance Refs. [21, 22, 23]. The analysis of Section and 



Section [2.3| has shown that there are two main factors producing differences between the line profiles 
emitted from the disc extending down within a boson star or the usually considered one accreting around 
a Schwarzschild black hole: 

1. The gravitational potential of a boson star is less deeper than that of a Schwarzschild black hole in 
the innermost region (r < 15r g ). 

2. The rotation velocity of the disc around a boson star is smaller than that of the disc around a 
Schwarzschild black at the inner region with r < 25r g 

Figures |3|j7| illustrate the line profiles in these two different situations using a completely relativistic 



calculation with a large range of parameters. Details of the code used can be found in references [17, 19]. 
These profiles correspond to the inner disc annulus with. T[ n = 6r g and r out = 30r g (Fig. || for a = 2 and 
Fig. |H for a = 3) or with r in = 2r g and r out = 30r 9 (Fig. || for a = 2 and Fig. for a = 3) in the case of a 
BS-AD (solid lines); or with rj n = 6r g and r out = 30r 9 in the case of a SchBH-AD (dotted lines). All the 
lines are normalized by their line flux. 



These are the main features of our results: 

• The lines emitted from the inner region of a BS-AD are somewhat narrower than those emitted from 
the inner region of a SchBH-AD. The red wing of the lines emitted from the inner region of a BS-AD 
does not extend to extremely low energy as that from a SchBH-AD because 1) the gravitational 
shifting of line photons (a factor of \/B(r)) emitted from the inner region of BS-AD is much less 
than that from the inner region of SchBH-AD (see Fig. [I] and 2) the Doppler redshifts of line photons 
in the case of BS-AD is also less than that in the case of SchBH-AD, since the rotating velocity of 
the inner disc is much smaller in the former case (see also Fig. |3|). The high energy blueward extent 
of the lines emitted from the inner region of BS-AD is determined by the inclination angle of the 
disc, which is similar to the case of SchBH-AD. However, for the same inclination angle, in the case 
of BS-AD, the blue horn extends to a larger energy than that in the case of SchBH-AD, since the 
former case suffers less gravitational redshifts. 
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• Generally, the line profiles emitted from BS-AD appears to have more peaks than that emitted from 
SchBH-AD (see Figs. ^Q). For very small inclination angles (e.g. 1°, see upper left panel in Figs. |5|j7]). 
the line emitted from the inner region of BS-AD is double peaked, and the low energy (more intense) 
peak is due to the emission from the inner edge where the gravitational shifting is the dominant 
factor (the transverse Doppler shifting is not large at the inner edge). The line from the inner region 
of SchBH-AD, on the contrary, tends to a single broad peak (and less intense) with skewed red wings 
due to much deeper gravitational potential and larger rotating velocity of the disc material near the 
marginal stable orbit. 

• In discs with large inclination angles (e.g. > 15°), the line profiles emitted from the inner region of 
BS-AD may have more peaks. A typical case with inclination angle O = 30° is illustrated in the left 
bottom panel in Fig. |4|: the profile for the inner region of a BS-AD has totally four local maxima 
compared to the only two local maxima for the case of SchBH-AD. From right to left, the first and 
the third maximum are due to the line emission in the outer part of the annulus [with r(v ~ w max )]- 
There are small displacements between these two maxima and the corresponding two maxima in the 
case of SchBH-AD. The second and the fourth maximum are due to the line emission from the inner 
part of the annulus [with r(v Zi v max )]. There are no such two corresponding peaks in the case of 
SchBH-AD because combined effects of the much deeper gravitational potential of the black hole and 
the larger disc rotation velocity shifts the line photons to lower or higher energy and smooth the line 
profile to the skewed one. 



4 Concluding remarks 

Observations of the central regions of active galaxies already exist. For some of them, the Iron Ko profile 
could be measured and interesting conclusions can be already drawn. For instance, given the observations of 



MCG-6-30-15, in which the red wing of the Fe Ka extends down to energy of 2-3 keV [24|, the computations 
of this paper have clearly shown that the central object in this galaxy cannot be a non-rotating massive 
boson star. It would be of interest, though, to study more complicated cases of accretion discs, particularly 
those around rotating boson stars. The observational data will be very much improved with forthcoming 
technologies. The NASA Constellation-X [^5|] mission, to be launched in 2008, is optimized to study the 
iron Ka line features and will determine the black hole mass and spin for a large number of systems. A 
detection of several peaks at lower energies, or of intense double peaks in faced on galaxies could be an 
impressive indication of a boson star nucleus. More importantly, this paper gives a new way to examine 
an interesting theoretical construct with real astrophysical data, helping to elevate the boson star concept 
to a falsifiable level. 
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Figure 4: Line profiles computed using complete relativistic method as described in the text. The line 
emitting region for all the cases is fixed to be the annulus with r; n = 6r g and r out 



30r g . The line 



emissivity law is set to be e oc r and the inclination angles are labeled in each panel. The solid lines 
and dashed lines represent the profiles of lines from boson star-accretion disc and from Schwarzschild black 
hole-accretion disc, respectively. 
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Figure 5: Same as Fig. |], but adopt a different emissivity law of e oc r 
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Figure 6: Same as Fig. ||, but adopt a different inner radius r in = 2r g of the line emitting annulus for the 
case of the boson star-accretion disc. 
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Figure 7: Same as Fig. B, but adopt a different emissivity law of e oc r 
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